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Density functional molecular dynamics simulations have been carried out to understand the finite 
temperature behavior of Auig and AU20 clusters. AU20 has been reported to be a unique molecule 
having tetrahedral geometry, a large HOMO-LUMO energy gap and an atomic packing similar to 
that of the bulk gold (J. Li et al., Science, 299 864, 2003). Our results show that the geometry 
of Auig is exactly identical to that of AU20 with one missing corner atom (called as vacancy). 
Surprisingly, our calculated heat capacities for this nearly identical pair of gold cluster exhibit 
dramatic differences. AU20 undergoes a clear and distinct solid like to liquid like transition with a 
sharp peak in the heat capacity curve around 770 K. On the other hand, Auig has a broad and flat 
heat capacity curve with continuous melting transition. This continuous melting transition turns 
out to be a consequence of a process involving series of atomic rearrangements along the surface 
to fill in the missing corner atom. This results in a restricted diffusive motion of atoms along the 
surface of Auig between 650 K to 900 K during which the shape of the ground state geometry is 
retained. In contrast, the tetrahedral structure of AU20 is destroyed around 800 K, and the cluster 
is clearly in a liquid like state above 1000 K. Thus, this work clearly demonstrates that (i) the gold 
clusters exhibit size sensitive variations in the heat capacity curves and (ii) the broad and continuous 
melting transition in a cluster, a feature which has so far been attributed to the disorder or absence 
of symmetry in the system, can also be a consequence of a defect (absence of a cap atom) in the 
structure. 



I. INTRODUCTION 

Recently discovered clusters and nanostructures of 
gold are found to have a rich chemistry with potential 
applications in materials science^ medicine^ and in the 
area of catalysis^ In particular, small clusters of gold 
have attracted interest as tips and contacts in molecular 
electronic circuits^ and also as chemical catalysts. 6 Ex- 
perimentally, even a small cluster such as Aus has been 
reported to catalyze the oxidation reaction of CO. 7 How- 
ever, these properties are reported to have strong size sen- 
sitive variations. Another factor influencing the applica- 
tion of gold clusters is their thermal stability. It is noted 
that several Au clusters undergo structural transforma- 
tion or tend to grow readily by migrating and merging 8 
under high-temperature conditions (500 K and above). 
These effects have important consequences in the appli- 
cations involving elevated temperatures and the growth 
mechanisms of clusters. In this context, a study on the 
finite temperature behavior of Au clusters is of consider- 
able importance. 

Since the pioneering reports on the possible applica- 
tions of gold clusters, a large amount of experimental 9 
and theoretical worksi - have been devoted to understand 
the structural and electronic properties of Au„ (n < 50) 
clusters. These reports have demonstrated that gold clus- 
ters have very different physical and chemical properties 
as compared to their bulk counterpart. A more recent 
exciting report has shown photo electron spectroscopic 



evidence of hollow golden cages with an average diameter 
of 5.5 A in the in the size range of 16 to 18 atoms^ These 
predictions were further supported by the theoretical cal- 
culations in the same report. However, AU20 is the most 
intriguing gold cluster reported so farJ^ Experimental 
studie s 11 ' 12 report this cluster to have a pyramidal struc- 
ture (tetrahedral symmetry) with each of the four faces 
representing the (111) surface of the Face Centered Cubic 
(FCC) gold. It is reported to have a large energy gap be- 
tween the Highest Occupied Molecular Orbital (HOMO) 
and Lowest Occupied Molecular Orbital (LUMO). This 
energy gap is greater than that of Ceo suggesting it to be 
highly stable and chemically inert. On the other hand, 
its structure with high surface area and large fraction of 
corner sites with low atomic coordination is expected to 
provide ideal surface sites to bind various molecules such 
as CO, O2 and CO2 for catalysis. The structure of AU19 
is also seen to be very similar to that of AU20 with one 
missing corner atom. In this context, it is interesting to 
have an understanding on the thermal stability of these 
two gold clusters having an atomic packing similar to that 
of bulk gold. Hence, in the present work, we study the 
finite temperature behavior of Auig and AU20 using the 
first principles Molecular Dynamics (MD) simulations. 

While several experimental and theoretical studies 
have been devoted to understand the ground state ge- 
ometries and chemical reactivity of gold clusters, there 
are very few reports on the finite temperature proper- 
ties of gold clusters ! 13 ' 14 The classical MD simulations by 
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Landman and co-worker s 13 ' 14 on medium sized Au clus- 
ters (150-1500 atoms) indicated that the clusters in this 
size range undergo a solid-to-solid structural transfor- 
mation around 700 K, before eventually melting around 
780 K. This is at a much lower value as compared to the 
bulk melting temperature of 1377 K. However, to the best 
of our knowledge, only one ab initio molecular dynamics 
study 15 attempting to understand the finite temperature 
behavior of Au clusters has been reported so far. 

As the cluster size reduces, the electronic effects 
play a more explicit role in controlling the struc- 
tural and thermal properties of the clusters. This is 
amply demonstrated by several first principles molec- 
ular dynamics simulationS ) 16 i 17 ' 18 i 19 ' 20 i 21 ' 22 ' 23 ' 24 which 
have successfully explained various experimental find- 
1^925,26,27,28,29,30,31,32,33 on the finite temperature be- 
havior of sodium, tin, gallium and aluminum clusters. 
These experimental studies have brought out various in- 
teresting phenomena such as higher than bulk melting 
temperatures in Ga and Sn clusters , 29 ' 30 and strong size 
dependent variations in the melting temperatures of Ga 
and Al clusters. 31 ' 32 However, the most surprising ex- 
perimental finding is the size sensitive behavior of the 
shape of the heat capacities where addition of even one 
atom is seen to result in a dramatic change of shape, 
prompting some of the clusters to be called as "Magic 
Melters" 4^ This means that while some clusters do un- 
dergo a conventional and clear melting transition, others 
undergo a near continuous transition making it very dif- 
ficult to identify any meaningful transition temperature. 
In a recent communication,— we have clearly demon- 
strated that a cluster with local "order" (an island of 
atoms connected with equal bond strengths) displays a 
well characterized melting transition with a distinct peak 
in the heat capacity curve, while a "disordered" cluster is 
seen to undergo a continuous transition with a flat heat 
capacity curve. Further, it is noted that this size sensi- 
tive nature in small clusters is related to the evolutionary 
pattern seen in their ground states and is seen to exist in 
clusters of sodium, gallium and aluminum. 24 ' 34 ' 35 

In what follows, we show that this dramatic variation 
in the shape of heat capacity curve is also observed in the 
present pair of gold clusters viz., Aui9 and Au2o- This 
observation has also thrown light on additional factors 
responsible for a continuous melting transition in clus- 
ters. As we shall see, in contrast to gallium or aluminum 
clusters, the flat or broad heat capacity curve in Auig 
is attributed to a "vacancy" in the surface (or a surface 
defect). This "vacancy" results in a chain of atomic re- 
arrangements leading to a restricted diffusion of atoms 
along the surface. In contrast, AU20 undergoes a rela- 
tively sharp melting transition with a clear peak in the 
heat capacity curve around 770 K. Thus, these results 
not only bring out an understanding on additional factors 
contributing to the broad melting transition in clusters 
but also on the relative thermal stability of these unique 
tetrahedral gold clusters. 



II. COMPUTATIONAL DETAILS 

We have optimized about 300 geometries for each of the 
cluster, to obtain the ground state geometry and several 
low energy isomers. The initial configurations for the op- 
timization were obtained by carrying out a constant tem- 
perature dynamics of 100 ps each at various temperatures 
between 400 to 1600 K. Once the ground state geometry 
is obtained, thermodynamic simulations are performed 
using Born-Oppcnhcimcr MD based on the Kohn-Sham 
formulation of Density Functional Theory (DFT)i 3 ^ The 
ionic phase space of the clusters is sampled classically in 
a canonical ensemble according to the method proposed 
by N6se<21 The MD simulations have been carried out 
using Vanderbilt's ultra soft pseudo potentials 3 ^ within 
the Local Density Approximation (LDA) for describing 
the core- valance interactions as implemented in the vasp 
package. 39 Energy cutoff of 13.21 Ry is used for the plane 
wave expansion of Au. We have used cubic super cells of 
length 20 A and have ensured that the results converge 
with respect to further increase in the energy cutoff and 
size of the simulation box. 

In order to have a reliable sampling, we split the total 
temperature range from 400-1600 K into at least 15 dif- 
ferent temperatures for both cluster sizes. We maintain 
the cluster at each temperature for a period of at least 
70 ps after equilibration, leading to a total simulation 
time of around 1 ns. Following the finite-temperature 
study, the ionic heat capacity of each cluster is computed 
using the multiple-histogram (MH) method . 40 ' 41 Vari- 
ous other thermodynamic indicators such as the mean- 
square displacements (MSD)'s of ions and the root-mean- 
square bond-length fluctuation (RMS-BLF or <5 rms ) are 
also computed. For the sake of completeness we briefly 
discuss these parameters. The parameter S rms is a mea- 
sure of the fluctuations in the bond lengths averaged over 
all the atoms and over the total time span. It is defined 
as 
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where N is the number of atoms in the system, r^- is the 
distance between atoms i and j, and (. . .)t denotes a time 
average over the entire trajectory. The MSD is another 
widely used parameter for analyzing a solid- like-to- liquid- 
like transition. In the present work, we calculate the 
mean square displacement's for individual atoms which 
is defined as 



<r?W> 



1 M 

■jj ^ [R-j(iom + t) — R/(iom)r 



(2) 



where Rj is the position of the Jth atom and we average 
over M different time origins to m spanning the entire tra- 
jectory. The MSD indicates the displacement of atom in 
the cluster as a function of time. In the solid-like region, 
all atoms perform oscillatory motion about fixed points 



(a) E=0.000 cV (b) E=0.065 eV (c) E=0.077 cV 




(d) E=0.113 eV (e) E=0.188 eV (f) E=0.212 eV 



FIG. 1: The ground state geometry and low lying isomers of 
AU19. The energy below is the difference in total energy of 
the isomer with respect to the ground state energy. 

resulting in a negligible MSD's of individual atoms from 
their equilibrium positions. In a liquid-like state, on the 
other hand, atoms diffuse throughout the cluster and the 
MSD's eventually reach a saturated value of the order of 
the square of the cluster radius. More technical details 
concerning the extraction of thermodynamic averages, in- 
dicators and computation of the heat capacity curve can 
be found in previous work4i 



III. RESULTS AND DISCUSSION 

We begin with a discussion on the ground state geome- 
tries and some representative low lying isomers of Auig 
and AU20 which are shown in the Fig. Q] and Fig. re- 
spectively. It is clearly seen from these figures that the 
ground state geometry of AU20 (Fig- I2H a )) i s a tetrahe- 
dron. The ground state geometry of Auig (Fig. [U-(a)) 
differs from that of AU20 by a single missing vertex atom 
of the tetrahedron. This is in agreement with the recent 
experimental and theoretical predictions. ^ Quite clearly 
both the structures are symmetric, with ordered trian- 
gles stacked over one another. The rest of the geometric 
parameters such as bond lengths, bond angles and dihe- 
dral angles are almost identical in both the ground state 
geometries. Thus, Auig can be considered as AU20 with 
a vertex defect. 

The atoms in Auig as well as AU20 are bonded to 
their first nearest neighbors with bond distances of ei- 
ther 2.63 A (shortest bonds in the cluster) or 2.75 A (next 
shortest bonds). It is interesting to display the connec- 
tivity of the shortest bonds to bring out the differences 
in them. Fig. 02-(a) an d Fig- I31~(b) show the distribu- 
tion of shortest bonds (2.63 A) for the case of Auig and 
Au2o, respectively. Clearly, the shortest bonds are dis- 
tributed only along the surface of both the clusters and 
form a closed network in AU20, while in Auig they form 



(a) E=0.000 cV (b) E=0.444 eV (c) E=0.445 eV 




(d) E=0.523 eV (e) E=0.553 eV (f) E=0.587 eV 



FIG. 2: The ground state geometry and low lying isomers of 
Au2o- The energy below is the difference in total energy of 
the isomer with respect to the ground state energy. 




FIG. 3: The distribution of shortest bonds (2.63 A) in Auig 
and Au2o- Rest of the inter atomic bond distances (for the 
first nearest nearest neighbor) in both the clusters are 2.75 A. 



an open network due to the missing atom. It turns out, 
that the presence of the vertex defect and open skeleton 
of shortest bonds in Auig play a significant role in the 
finite temperature behavior of the cluster around 650 K 
and initiate a set of restricted atomic rearrangements on 
the surface. 

The low lying isomers of Auig (shown in Fig. [l]-(b) 
to Fig. [I]-(f)) are clearly devoid of a regular triangular 
arrangement of atoms seen in the ground state configu- 
ration. The first low lying isomer of Auig (Fig. [T]-(b)) 
is nearly 0.065 eV higher in energy as compared to the 
ground state configuration. Auig has several isomers 
with continuous energy distribution between 0.065 eV 
to 0.7 eV some of which are shown in Fig. UK ) to Fig. 
[T]-(f). This can be contrasted with first low lying isomer 
of AU20 (shown in Fig. HJ-(b)) which is almost 0.44 eV 
higher than the ground state geometry. This structure 
has one central atom and rest of the 19 atoms arrange 
around this central atom so as to have a highly deformed 
tetrahedron. This configuration is degenerate with the 
hollow cage configuration of AU20 (Fig. H]-(c)). The fact 
that there are a couple of isomers well separated from the 
ground state geometry correlates very well with the ex- 
istence of a relatively sharp and well defined peak in the 
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FIG. 4: The Heat-Capacity Curves of Auig and Au2o- 




(d) (e) (f) 

FIG. 5: (a-e) Snapshots of restricted rearrangement of atoms 
in Auig around 650 K. (f) Arrows depict the continuous 
atomic rearrangements that take place to fill in the missing 
cap atom between 700 K-900 K. 



heat capacity curve of Au^o- Some other representative 
low lying isomers of AU20 are shown in Fig.[2]-(d) to Fig. 
[2Hf). We also note that the atoms in all the high energy 
configurations of both clusters are bonded to each other 
through a much wider and continuous range of Au-Au 
bond lengths (ranging between 2.67 A-2.90 A) as com- 
pared to those in the ground state geometry. 

Now, we present the finite temperature behavior of 



both the clusters, we begin with a discussion on the 
calculated heat capacity curves which is shown in Fig. 2] 
The figure brings out a remarkable feature, viz., a signif- 
icant size sensitivity nature of the heat capacity curves. 
The heat capacity curve of AU20 has a clear and recogniz- 
able peak around 770 K, with a width of about 250 K. 
In contrast, the heat-capacity curve of Auig, a cluster 
with a vacancy, shows a broad and almost continuous 
solid-to-liquid transition between 650 K-1200 K. Thus, 
this is yet an another example of dramatic change in the 
shape of the heat capacity curve with the addition of a 
single atom. As already noted this size sensitive nature 
has been observed earlier in Ga, Al clusters experimen- 
tallyM 

It is possible to make a detailed analysis of the ionic 
motion by examining the trajectories of the clusters. An 
analysis of the ionic motions of Auig reveals the cluster 
to vibrate around its ground state geometry until 600 K. 
Around 650 K, the cluster undergoes a peculiar struc- 
tural rearrangement so as to fill in the vacancy (the apex 
atom). The snapshots of this structural rearrangement 
are shown from Fig.[5]-(a) to Fig.[5]-(e). Fig.[5]-(a) shows 
Auig with a missing cap atom. Note that, at the end 
of the structural rearrangement (Fig.[5]-(e)), the vacancy 
which is present on the top in Fig. El-(a) is shifted to the 
bottom edge. We denote the edge consisting of atoms 
'A', 'B' and 'C as the reference edge. Coming to the 
details of the structural rearrangement, as the cluster 
evolves around 650 K, it is seen that these atoms in the 
reference edge push themselves upward (see Fig. [5]-(b) 
to Fig. [S]-(d)). The rest of atoms in the cluster undergo 
minor displacements around their equilibrium positions 
during this process. At the end of this rearrangement 
(see Fig. [5]-(e)), it is seen that the top edge atom (atom 
'C') in Auig moves to cap the missing vertex atom seen 
in Fig. [5]-(a). The next edge atom ('B') moves up to 
occupy the position initially occupied by atom 'C and 
atom the 'A' occupies the position occupied by 'B'. Now 
this creates a vacancy or defect at the position initially 
occupied by the atom 'A'. Thus, we now have a Auig 
cluster which is rotated by 90 degrees in the anti clock 
wise direction (Fig. 0-(e)) with respect to Fig. [S]-(a). 

Around 650 K, only a single edge is displaced so as 
to cap the missing vertex atom. Between 700 K-900 K, 
we see a continuous displacement of atoms along all the 
edges as shown in Fig. [S]-(f) as the vacancy is shifted 
from one vertex to the other vertex. A remarkable fea- 
ture of this motion is that the overall shape of the cluster 
remains approximately tetrahedron with a missing cap. 
Around 1000 K, the tetrahedron structure is destroyed 
and the cluster visits its first and second high energy con- 
figurations. The cluster finally melts completely above 
1200 K. This leads to a broad feature (between 650 K- 
1200 K) in its heat capacity curve. 

In contrast, the ionic motion of AU20 shows all the 
atoms to vibrate around their initial positions until about 
750 K. The cluster undergoes a structural transformation 
from the ground state geometry to the first isomer shown 
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in Fig. [2]-(b) around 800 K. The cluster visits other iso- 
mers around 900 K and melts completely above 1000 K 
leading to a clear and relatively narrow melting transi- 
tion. 

This contrasting behavior is brought out more clearly 
by examining the MSD's of the individual atoms. In 
Fig. [6l we show the MSD's of individual atoms in both 
clusters in the temperature range of 450 K-1000 K. It is 
clearly from Fig.[|)]-(a) and Fig.[5]-(b) that atoms in both 
the clusters vibrate around their equilibrium positions 
at 450 K. Around 650 K, the rearrangement of atoms 
along one edge in Au 19 is reflected in slightly higher mean 
square displacements (around 3 A) in Fig. [6]-(c). These 
values increase continuously in Auig as shown by the 
typical behavior around 800 K. In contrast, the MSD 
values in AU20 are negligible until 750 K (< 0.5 A) and 
the values increase sharply around 800 K. It is precisely 
at this temperature that the tetrahedron gets destroyed. 
The MSD's of Auig and AU20 saturate around 35 Aat 
1200 K and 1000 K respectively indicating the presence 
of a liquid like state. 

This contrasting behavior is some what weakly re- 
flected in in average root-mean-square bond-length fluc- 
tuation (i rms ) of Auig and AU20, which is in any case a 
quantity averaged out over all the atoms. In Fig. [7| we 
show the <5 rms for both the clusters. As expected AU20 
shows a sharp transition indicated by a jump in 5 Tms value 
from 0.07 to 0.20. 



IV. SUMMARY AND CONCLUSIONS 





In this work, we have presented the results of first prin- 
ciple molecular dynamics simulations on Auig and AU20, 
clusters which have atomic packing similar to that of the 
bulk gold. In spite of the fact, that the geometry of both 
the clusters is nearly identical except for a single vertex 
atom, they exhibit dramatic differences in the shape of 
their heat capacity curves. We have shown that these 
differences are induced by vacancy. The vacancy in Auig 
induces a restricted diffusive motion along the surface 
of the cluster leading to a continuous melting transition. 
In contrast, AU20 exhibits a sudden and clear melting 
transition. It may be noted that such a size sensitive 
nature of the heat capacity curves has been observed ex- 
perimentally in Al and Ga clusters and in Na clusters 
during the ab initio molecular dynamic simulations. In 
these studies the size sensitive nature of the heat capac- 
ity curves was attributed to the nature of the "disorder" 
in the ground state geometry. The present work clearly 
shows that this size sensitive behavior is also driven by 
the vacancy in the other wise perfect and symmetric "lat- 
tice" . The work also demonstrates that the size sensitive 
variations in the melting characteristics to be generic in 
nature. Finally, the contrasting finite temperature be- 
havior reported in the present gold clusters could have 
several implications in the applications of these clusters 
and is a topic of further research interest. 



FIG. 6: The Root Mean Square Displacements of atoms with 
respect to the simulation time (ps) in Auig at various tem- 
peratures. 
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